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ABSTRACT 

We present the first results of an analysis of the properties of the molecular gas in the 
nuclear regions (r < 300 pc) of a sample of six nearby galaxies, based on new high spatial res- 
olution observations obtained in the /iT-band with the near-infrared integral field spectrograph 
SINFONI at the Very Large Telescope. We derive two-dimensional distributions of the warm 
molecular and ionized gas from the H2, Bry and Hei emission lines present in the spectra of 
the galaxies. We find a range of morphologies, including bar- and ring-like distributions and 
either centrally peaked or off-centre emission. The morphologies of the molecular and the 
ionized gas are not necessarily coincident. The observed emission-line ratios point towards 
thermal processes as the principal mechanism responsible for the H2 excitation in the nuclear 
and circumnuclear regions of the galaxies, independently of the presence of an active nucleus. 
We find that a rescaling of the H2 2.12 fim emission-line luminosity by a factor fi - 1200 gives 
a good estimate (within a factor of 2) of the total (cold) molecular gas mass. The galaxies of 
the sample contain large quantities of molecular gas in their centres, with total masses in the 
~ 10 5 - 10 s M0 range. Never the less, these masses correspond to less than 3 per cent of 
the stellar masses derived for the galaxies in these regions, indicating that the presence of gas 
should not affect black hole mass estimates based on the dynamical modelling of the stars. The 
high-spatial resolution provided by the SINFONI data allowed us to resolve a circumnuclear 
ring (with a radius of ~ 270 pc) in the galaxy NGC 4536. The measured values of the Bry 
equivalent width and the He i/Bry emission-line ratio suggests that bursts of star formation 
occurred throughout this ring as recently as 6.5 Myr ago. 

Key words: galaxies: nuclei - infrared: galaxies - galaxies: ISM - ISM: molecules - line: 
formation 



1 INTRODUCTION 

It is widely accepted that supermassive black holes (SMBHs), 
with masses between 10 6 - 10 10 Mq, sit at the centre of el- 
liptical galaxies and bulges of spiral galaxies. The observational 
correlation between the mass of the central black hole (M BH ) 
and the bulge mass or the luminosity and the velocity disper- 
sion (<r) of the galaxy's bulge, strongly suggests a connection be- 
tween the formation of the central black holes and the formation 



1 1 1 n 

and evolution of their hosts (e.g., Kormendv & Richstone 1995 


Maaorrian et al. 1998 


; Ferrarese & Merritt 2000; Gebhardt et al 


2000t ITremaine et al .1 


20021; IMarconi & Huntll2003l; lHarine & Rix 
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l2004lGultekin et al.ll2009T) . However, the universality of these re- 
lationships has not yet been firmly established. This is mainly due 
to the fact that they are still based on a limited sample of galaxies, 
biased towards early types. 

In this spirit, an observational spectroscopic survey of nearby 
galaxies, u sing the near-infrared (NIR) integral field s pectrograph 
SINFONI jEisenhauer et alj|2003l : bonnet et alj|2004h at the Very 
Large Telescope (VLT), was undertaken. The survey is aimed at in- 
creasing the number of M B h estimates to explore the M BH -bulge 
relations, especially at the poorly known low- and high-mass ends, 
and testing the BH-bulge formation scenarios. This paper is part of 
a series reporting the results of the SINFONI program (Nowak et 
al. 2007, 2008, 2010; Rusli et al. 2011, 2012a, 2012b; Bender et al. 
in preparation; Erwin et al. in preparation; Saglia et al. in prepara- 



© 0000 RAS 



2 Mazzalay et al. 



tion). A thorough analysis and derivation of black-hole masses via 
extensive stellar dynamical modelling can be found in these papers. 

Of the total sample of galaxies observed with SINFONI, about 
one fourth exhibit signatures of gas emission in their nuclei. The 
gaseous material is considered to be the primary fuel source of 
the nuclear activity. It is not only necessary for the formation and 
growth of the nuclear BHs but it is also a fundamental ingredient for 
the nuclear and circumnuclear starburst activity of galaxies. How is 
the large amount of gas required to fuel the nuclear activity trans- 
ported to the inner regions of the galaxies? What is the interplay 
between star formation and central SMBHs? These are among the 
important questions that need to be answered in order to improve 
our understanding of galactic nuclei, and are currently the focus of 
intensive research. 

The fuel necessary to drive the nuclear activity has to be trans- 
ported from galactic scales (~ 10 kpc) to scales of less than 1 kpc 
and into the nucleus. Depending on the scales involved, a num- 
ber of mechanisms that remove angular momentum from the gas, 
thus bri nging it clos er to the centre, have been proposed (see re- 
view bv |jogeell2006l . and references therein). Large-scale bars (and 
the star-forming rings sometimes associated with them), as well 
as tidal interaction and mergers, seem to play a key role in the 
transport of the g as into the inner kpc. On smaller scales, nuclear 
nested bars (e.g., IShlosman et alj|l989t) and nucle ar spirals (e.g., 
lEnglmaier & Shlosmanll2000l : lMacieiewskill2004bllal) are some of 
the means of removing angular momentum that have been pro- 
posed. 

The study of the distribution and properties of the molecular 
gas and star-forming regions in galactic nuclei can provide clues on 
the mechanisms channelling the gas to the centre and thus fuelling 
nuclear activity, and can also help establish the link between the 
starburst and nuclear SMBHs. The SINFONI data set is ideal for 
this task, since it simultaneously provides high-spatial resolution 
and the advantages of NIR wavelengths, like low dust extinction 
and a number of emission lines that can be used as diagnostic for 
the properties of the molecular gas and the nuclear activity. 

This is the first of two papers presenting new integral field 
spectroscopic (IFS) data for 6 of the galaxies that show molecular 
emission in our SINFONI sample of 33 galaxiefl Both are aimed 
at characterizing the properties of the nuclear gas and determining 
its effect on (and possible use in) the estimation of the mass of the 
central BH. Here we focus on the emission line gas distribution and 
mass content of the galaxies, considering the implications for the 
estimation of BH masses. In a second paper (Mazzalay et al. 2012, 
in preparation) the kinematics of the gas is analysed and compared 
with that of the stars. The layout of this paper is as follows. In Sec- 
tion|2]we describe the new SINFONI observations and the data re- 
duction process. In Section[3]we describe the main properties of the 
galaxies in the sample and analyse the morphology of the nuclear 
regions. Section|4]presents the results on the physical conditions of 
the molecular and ionized gas. An inventory of the total (gas plus 
star) mass content is given in Section [5] A summary and our main 
conclusions can be found in Section[6] 



1 Two further galaxies, N GC 1316 and NGC 3368, s how molecular emis- 
sion and were discussed bv lNowak et all <2008ll20ld) . 



2 OBSERVATIONS AND ANALYSIS OF THE DATA 
2.1 IFS observations 

The data presented here were gathered as part of our observational 
campaign to measure SMBH masses using the NIR integral-field 
spectrograph SINFONI at the 8 m Very Large Telescope UT4. The 
galaxies were observed in several runs between 2007 and 2009. A 
short description of the observations can be found in Table [T] We 
covered the 1.95-2.45 yum wavelength range with the A"-band grat- 
ing using two different spatial samplings: a high-spatial resolution 
scale of 0.05 x 0.10 arcsec pixeL 1 (hereafter HR data) and a low- 
resolution scale of 0. 125x0.250 arcsec 2 pixeL 1 (hereafter LR data). 
The associated fields of view (FOVs) are approximately 3" x3" and 
8" x 8", respectively. The observations were performed following 
a standard object-sky-object strategy, with individual exposures of 
10 min dithered by a few spatial pixels (or 'spaxels'). This pro- 
cedure allows a good correction of bad pixels and cosmic ray re- 
moval. The total on-source exposure times are listed in column 4 of 
TableQ] The high-spatial resolution observations were AO-assisted, 
emplo ying either the laser guide star (LGS) PARSEC (Rab ien et alj 
|2004) or the natural guide star (NGS) mode using the galaxy nu- 
cleus, while the low-resolution observations were taken in natural 
seeing. In order to assess the AO performance and characterize the 
point spread function (PSF) due to atmospheric turbulence, we reg- 
ularly observed a PSF star immediately after the science exposure 
sequence, in addition to the normal (non-AO) observation of a tel- 
luric standard star. The PSF stars for individual galaxy observations 
were chosen to have optical (B - R) colours and K-band magni- 
tudes similar to that of the inner r < 3" of the galaxy nuclei, in 
order to achieve a similar AO correction. Column 6 of TableQ]lists 
the FWHM of the PSF star associated with each AO-assisted ob- 
servation and, for the cases of the non-AO observations, the seeing 
values at 2.2 jim estimated from those given by the Differential Im- 
age Motion Monitor (DIMM). Note that, since the DIMM seeing is 
measured at 0.5 jim, we assumed that the resolution 6 depends on 
the wavelength as 9(A) ~ A~ 1/5 in order to have an estimation of the 
seeing at 2.2 jim. Both high- and low-spatial resolution configura- 
tions give a spectral resolution R ~ 4000 at 2. 12 /urn, corresponding 
to an instrumental dispersion of cr ins ~ 30 km s . 



2.2 Data reduction 

The dat a were reduced using a custom pipel ine incorporating ES- 
OREX dModigliani et al.l2007h and SPRED dSchreiber et alj2004 
lAbuter et al.ll2006l) . The reduction process comprised all standard 
reduction steps, including bias subtraction, flat-fielding, bad pixel 
removal, detector disto rtion and wavelength calibrations, sky sub- 
traction jDaviesI |2007h . reconstruction of the object data cubes 
and, finally, telluric and flux calibrations. Through the reduction 
tasks, the data cubes were re-sampled to a spatial scale of 0.05 x 
0.05 arcsec 2 pixeL 1 and 0.125 x 0.125 arcsec 2 pixeL 1 for the HR 
and LR data, respectively. Telluric correction and flux calibration 
of the reconstructed data cubes were performed using the telluric 
standard stars observed at the beginning or the end of the corre- 
sponding observing block at an airmass similar to that of the galaxy. 
Each star spectrum was extracted from an aperture with a diame- 
ter of 5-7 times the mean FWHM and its continuum was modelled 
with a blackbody curve of the same effective temperature and K 
magnitude as the star. To test our flux calibration, we used 2MASS 
images in combination with high resolution HST (NIC2 and NIC3) 
and William Herschel Telescope (WHT) NIR images available for 



© 0000 RAS, MNRAS 000, 000-000 



Molecular gas in the centre of galaxies - I 3 



Table 1. Log of observations. The columns show the field of view (FOV) covered by the observations, observing date, on-source exposure time, AO mode 
(where NGS stands for 'natural guide star' and LGS for 'laser guide star' ), estimated FWHM of the point spread function and position angle of the observation. 



Galaxy 


FOV 


Obs. date 


T exp [min] 


AO mode 


PSF 


Obs. PA. 


NGC 3351 


3" X 3" 


2009-04-19 


40 


LGS 


0.20" 


10 




3" x 3" 


2009-04-20 


40 


LGS 


0.18" 


10 




8" x 8" 


2009-04-21 


40 


no AO 


0.45" 


10 


NGC 3627 


3" x 3" 


2007-03-21 


80 


LGS 


0.24" 


-12 




3" x 3" 


2007-03-25 


50 


LGS 


0.27" 


-12 




8" x 8" (N) 


2007-03-25 


80 


no AO 


0.67" 


-12 




8" x 8" (S) 


2009-04-22 


20 


no AO 


0.67" 


-12 




8" x 8" (S) 


2009-05-18 


20 


no AO 


0.59" 


-12 


NGC 4501 


3" x 3" 


2008-03-12 


80 


NGS 


0.13" 





NGC 4536 


3" x 3" 


2009-04-19 


60 


LGS 


0.18" 


125 




3" x 3" 


2009-04-22 


80 


LGS 


0.18" 


125 




8" x 8" 


2009-04-21 


80 


no AO 


0.74" 


125 


NGC 4569 


3" x 3" 


2008-03-08 


80 


NGS 


0.15" 


25 


NGC 4579 


3" x 3" 


2008-03-08 


40 


NGS 


0.15" 


96 




3" x 3" 


2008-03-09 


80 


NGS 


0.15" 


96 



NGC 3351, NGC 4569 and NGC 4579. The images were obtained 
from the HST archi ve, NASA/IPAC Infrared Science archive and 
iKnapen et alj d2003l) . We find an agreement between the different 
data sets consistent to 20 per cent. 

2.3 Construction of the 2D maps 

The SINFONI data presented here allow us to study the morphol- 
ogy and physical conditions of the circumnuclear molecular gas 
of the galaxies in the sample. To improve the signal-to-noise ratio 
(SNR) of the two-dimensional data and, at the same time, preserve 
the highest possib le spatial resolution, we u sed the adaptive 2D- 
binning method o f lCappellari & CopirJ I2003T) . This method makes 
use of Voronoi tessellations to bin pixels together by accreting new 
pixels into the bin until a desired SNR is reached. Here the SNR 
was defined as the ratio between the flux of the particular emission 
line analysed and the rms of the continuum, op t imally weighted as 
described in section 2. 1 o f lCappellari & Copml 120031) . The chosen 
SNR threshold depended on each individual galaxy and was se- 
lected so as not to compromise spatial resolution in the high SNR 
regions and, at the same time, be able to extract information from 
the regions which showed signs of emission but too low SNR to 
obtain otherwise meaningful line parameters. 

In order to extract information about the molecular gas we fit- 
ted a Gaussian component to the H2 2.12 /im emission line in the 
continuum-subtracted spectra of the galaxies, obtaining the central 
wavelength, width, and intensity of the line at each spatial bin. This 
line was selected since it is one of the strongest Hi lines observed in 
the K-bmd range and it is located in a region not contaminated by 
other spectral features. We used the measured intensity to construct 
emission-line flux distribution maps of the H2 emission line gas. 
The stellar continuum was fitted an d subtracted applying the penal - 
ized pixel fitting method (pPXF) of lCappellari & Emselleml d2004t) . 
For this we used a stellar template constructed as a positive linear 
combination of the spectra of six late-type stars observed with the 
same configuration as the galaxy data (in this way we ensure that 
the stellar template and the objects spectra have the same instru- 
mental broadening). Particular care was taken in masking regions 
contaminated by emission lines and/or spurious features (i.e., bad 
pixels or bad sky-lines subtraction). 

Visual inspection of the emission-line profiles showed no sig- 



nificant deviation from a single-Gaussian component for all the 
galaxies except NGC 4579. The H 2 emission-line profiles of this 
galaxy display strong asymmetries and double components in some 
regions. In consequence, to obtain a good estimation of the total 
integrated line-flux for this galaxy, we described its line profiles 
as the sum of two Gaussian components. The resulting maps of 
Hi 2.12 yum flux distribution are shown and analysed in the next 
section. 

Additionally, Bry emission was observed in the spectra of 
NGC 3351, NGC 3627 and NGC 4536, and the He 1 2.06 fim emis- 
sion line was present in the spectra of NGC 4536. Following the 
same procedure as in the case of H2, we constructed the flux distri- 
bution maps for these lines and included their analysis in the next 
section. 

As an example, Fig. [T] shows the integrated continuum- 
subtracted spectrum of the inner ~ 8" x 8" of NGC 4536 obtained 
in the LR configuration. 



3 MORPHOLOGY OF THE CIRCUMNUCLEAR GAS 

All the galaxies discussed in this paper display Hi emission lines 
in their spectra and, in some cases, Bry and He 1 2.06 /jm emission 
were also observed. In this section we describe and analyse the 
emission-line flux distribution maps of each individual galaxy, as 
well as relevant previous observations which will allow us to under- 
stand the context in which these data reside. All the SINFONI maps 
are centred at the A"-band continuum maximum; the continuum was 
estimated by integrating the signal in the 2.1-2.3 /im wavelength 
range. The main properties of the galaxies are summarised in Ta- 
ble^ The inclination and position angles of the galaxies were de- 
rived from a combination of analysis of the shape of outer galaxy 
isophotes and (where available) large-scale gas kinematics in the 
literature; details for individual galaxies will be presented in Erwin 
et al. (in preparation). 

3.1 NGC 3351 

NGC 3351 (M95) is one of the main members of the Leo I 
group of galaxies and it is classified as an early type barred 
spiral galaxy, with no signs of activity in its nucleus. It hosts 
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Figure 1. Continuum-subtracted SINFONI spectrum of the inner i 



' of NGC 4536. The labels mark the position of the emission lines identified. 



Table 2. Properties of the galaxies of the sample. 





NGC 3351 


NGC 3627 


NGC 4501 


NGC 4536 


NGC 4569 


NGC 4579 


Morphology " 


SB(r)b 


SAB(s)b 


SA(rs)b 


SAB(rs)bc 


SAB(rs)ab 


SAB(rs)b 


Nuclear activity 




LINER/Sy2 " 


Sy2" 


LLAGN b 


LLAGN c 


LINER/Syl.9 


Distance (Mpc) 


10.0 d 


10.1 d 


16.5 ' 


14.9 d 


16.5 e 


16.5 ' 


Linear scale (pc arcsec -1 ) 


49 


49 


80 


72 


80 


SO 


Inclination (deg) 


46 


65 c 


60 


69 


69 


40 


PA (deg) 


10 


178 / 


140 


125 


25 


95 


PA of the stellar bar (deg) 


112* 


161 '' 






15 * 


58; 



(n) Taken from NED. {b 



2007]). ( ^Casasola etal 



McAlpine et alj 
' ' 201 lh . ( * lErwir 



201 ll). (c ]Ho et all fl997l). '^Freedman et all feOOll) . M Mean Virgo cluster dist ance iMei et alj 
2005). ( * ]Sheth et aljj2002h . (i ljogee et al.lbOo'sli ^ lGarcfa-Burillo et alJfeOOSh . 



a star-forming nuclear ring with a major axis of ~15", which 
dominates the appearance of the inner region of this galaxy, 
and a larger ring of Hn regions that encircles its stellar bar 



Rubin et al. 



lg ot L in regions tha t encircles its stell ar bar 
j[l975l ; |Pevereuxetai]|l992l ; IColina et alj|l997l : 



iKnapen et alj 2002 ). Aligned almost perpendicular to the stellar 
ba r, a nuclear 'mo l ecular gas bar' ~ 19" long was first reported 
by iDevereux et al. I l ll992l) . Higher spatial resolution observations 
resolved this into two peaks of CO emission at a radius of ~ 7", 
coincident with t he circumnuclear ring observed in Ho- images 
jjogee et alhoOSl) . 

NGC 3351 was observed with SINFONI in two different con- 
figurations, with high- and low-resolution spatial samplings, cover- 
ing ~ 3" x 3" and ~ 8" x 8" FOVs. The H 2 and Bry flux distri- 
bution maps derived from these data sets are presented in Figure|2] 
The HR data shows that the nuclear emission is resolved into two 
different clumps embedded in diffuse emission (upper-left panel of 
Fig- EJ- The innermost emission is elongated along the major axis 
of the galaxy, with a relatively constant flux over the inner ~ 1" 
and a sharp drop to values below 20 per cent of the flux maximum 
further out. A smaller off-centre source is located at ~ 1" towards 
the NW of the nucleus. These two clumps are not individually re- 
solved in the LR data (upper-right panel of Fig. f2j, where we can 
see an almost symmetric central source connected towards the SW 
by a low-luminosity ridge to the inner edge of the nuclear ring ob- 
serve d in Ha and UV images dColina et alJI 19971 : IPlanesas et alj 
1997) of this galaxy. Although only part of the ring falls within the 
SINFONI field, the observed H 2 and Bry distribution found here 
agrees very well with previous observations of Ha emission-line 
maps, especially at the location of the hot spots. The brightest Bry 



emission can be asso ciated with the Ha hot spot labelled as R2 by 
IPlanesas et al. and the emission observed at the east edge of 

the SINFONI field can b e associated with the R5, R6 and R7 Ha 
hot spots (see figure 8 of lPlanesas etai]|l997t) . However, the Bry 
map does not show the nuclear emission observed in both the H 2 
maps presented here and previous Ha images. This is probably due 
to the intrinsic weakness of the Bry line, which falls below the sen- 
sitivity limit of our data. Hints of nuclear Bry emission is observed 
in the inner ~ 1" of the HR spectra, although it is not strong enough 
to allow the construction of an emission-line map. 

3.2 NGC 3627 

NGC 3627 (M66) is a high-inclination strongly barred galaxy. It 
belongs to the Leo Triplet group of galaxies, showing some signa- 
tures of past interaction with its neighbour galaxy NGC 3628. Op- 
tical emission-line ratios are inconclusive with respect to the type 
of nucleus this galaxy harbours. Whether it is a L INER/H II com- 
posite or Seyfert 2 like nucleus l lHo etalj[l997h . the presence of 
a variable, flat-spectrum nuclea r radio source is consistent with a 
low-luminosity AGN (LLAGN. IFilho et alj|200d) . Broad-band op- 
tical images of NGC 3627 reveal a pronounced spiral pattern and a 
central region dominated by several strai ght dust lanes that ru n par- 
allel to the host galaxy semi-major axis jMartini et al]|2003l) . The 
CO emission shows a nuclear peak and extends along the leading 
edges of the stellar bar, forming two broad p eaks at the bar en ds, 
from where the spiral arms trail off from (e.g. JSheth et alj2002h . 

Molecular and recombination lines of hydrogen were ob- 
served in the SINFONI spectra of NGC 3627. The left panels of 
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Figure 2. Flux distribution of H2 2.12 fim (upper panels) and Bry (lower panel) emission lines observed in NGC 3351. The panel on the left corresponds to the 
high-resolution data (FOV ~ 3" X 3") and the panels on the right correspond to the low-resolution data (FOV ~ 8" X 8")- The spatial orientation, indicated in 
the upper-left panel, is the same for all the panels. The position of the major axis of the galaxy (dashed line) and the orientation of the stellar bar (dotted line) 
are also indicated. 



Fig.[3]show the flux distribution maps of the H2 2.12 /im and Bry 
lines in the inner ~ 3" x 3" (HR maps) of this galaxy. Addition- 
ally, two sets of low-resolution observations were taken with a spa- 
tial offset along the -12° position angle, yielding a total FOV of 
~ 8" x 15". The final H2 distribution map obtained from these data 
(LR map) is presented in the right panel of Fig. [3] No evidence of 
Bry emission was found at distances larger than a few arcseconds 
from the nucleus, therefore we only show its distribution derived 
from the HR data. 

As it can be seen from Fig. [3] while the emission of Bry ap- 
pears strictly confined to the nucleus, within a region of < 50 pc in 
radius, the morphology of the H 2 emission in the inner 3" x3" is ex- 
tended along N-S direction. Further out, the gas is mainly located 
along a bar-like structure, with two hot-spots located approximately 
at distances of 6.5" and 4.5" from the nucleus along a position an- 



gle of ~ 17° (from north through east). Additionally, weak emission 
extending southward is observed beyond the southern hot-spot. The 
morphology displayed by the H2 emission line is r emarkably con- 
sistent with the 12 CO(2-l) morphology reported bv lCasasola et alj 
fcoill) for the inner 20" x 20" of NGC 3627 (see their figure 4). 
The position angle of the molecular bar-like structure observed in 
the innermost regions of NGC 3627 differs from the position angle 
of the stellar bar (PA = 161°, indicated with a d otted line in Fig.Q 
and the large-scale molecular bar (PA = 164°. ISheth et alj|2002h . 
This is because we are looking well within the transition region, 
where the molecular gas goes from leading the northern part of the 
stellar bar to lead the southern one (see, for example, figure 1 of 
lReganetal.ll200ll) . 
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Figure 3. Flux distribution maps of the H 2 and Bry emission lines observed in NGC 3627. The panels on the left correspond to the high-resolution data and the 
panel on the right corresponds to the low-resolution data. The spatial orientation, indicated in the upper-left panel, is the same for all the panels. The dashed 
line in the H2 maps indicates the position of the major axis of the galaxy and the dotted line corresponds to the stellar bar position angle. 



3.3 NGC 4501 

The Seyfert 2 galaxy NGC 4501 (M88) is one of the largest spi- 
ral gakixiej>JriJhe_Vh^ 6-band images presented 
bv lElmegreen & Elmegreer] i 19871) showed multiple external arms 
and two symmetric inner arms. In th e A"' -band, the spiral s truc- 
ture is continuous w ithin r < 3 kpc ( lElmegreenetal.lll999h and 
ICarollo et alj l l 19981) were able to observe spiral-like dust lanes 
(thought to indicate the location of shocks) down to the nucleus, 
using WFPC2 F606W images. High-resolution interferometric ob- 
servations of the 12 CO(1-0) emiss ion in the central 5 k pc of this 
galaxy were analysed in detail by lOnodera et alj J2004h . The CO 
maps showed spiral arms, associated with the dust lanes found by 
HST, and a strong central condensation of radius ~ 5" resolved 
into double peaks separated by 4.7", which are located on the 'root 
of the dust spirals'. Even though NGC 4501 does not show evi- 



dence for a bar, the gas concentration in its centre i s relatively high, 
and m ore typical of barred than unbarred galaxies l lSakamoto et alj 
Il999al) . 

Figure [4] shows the flux distribution map of the H2 2.12 /jm 
emission-line in the inner 3" x 3" of NGC 4501 observed with 
SINFONI. Two distinct structures can be resolved in the integrated 
flux map: an asymmetric nuclear component (note that the max- 
imum of the H 2 emission is offset from the continuum maximum 
indicated by the plus symbol) surrounded by two arcs that resemble 
an incomplete ring-like structure. The arcs traced by H 2 towards the 
NW and SE of the nucleus see ms to be located on the inside edges 
of the two peaks seen in CO bv lOnodera etall J2004I) . indicating a 
gradient in the temperature of the molecular gas in the sense that the 
warm molecular gas (traced by the H 2 2. 12 yum emission) is located 
closer to the nucleus than the cold molecular gas (traced by the 
CO emission). A comparison between the H 2 emitting-gas distri- 
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Figure 4. Flux distribution map of the H2 emission line observed in 
NGC 4501. The dashed line represents the orientation of the major axis 
of the galaxy. 
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Figure 5. Left: Unsharp mask of the HST/WFPC2 F547M image of the 
inner 10" X 10" of NGC 4501. Dust lanes trace spiral patterns around and 
to the nucleus. Right: Zoom into the inner 3" X 3" with the H2 contours 
overlaid. 



button and the nuclear spiral-like dust lanes exhibit by NGC 4501 
can be seen in Fig [5] This figure shows an unsharp mask of the 
HST/WFPC2 F547M image (retrieved from the HST archive) of 
the inner regions of the galaxy, with the H 2 contours overlaid. As 
in the case of the CO emission, a close spatial coincidence between 
the H 2 molecular gas and the dust lanes is observed, especially in 
the north region where both the molecular and dust structures are 
better defined, suggesting a close relation between the two media 
down to distances of less than 1" from the nucleus. 



3.4 NGC 4536 

NGC 4536 is a large spiral galaxy located in the Virgo clus- 
ter. While it has been classified as a SA B(rs)bc galaxy (e.g., 
Ide Vaucouleurs et"al]|l99ll ; iButa et a! the evidence for the 

presence of a bar is not yet clear (see for example Pompea & Rieke 
1990; Shaw et al. 1995, Erwin et al. in preparation). Alth ough op- 
tical l ine ratios place this galaxy in the H 11 r egion range jHo et alj 
1 1997b . the presence of hi gh-ionization lines dSatvapal et alj|2008h 
and its X-ray spectrum ( IMcAlpine et al J 1201 lh suggest that this 



galaxy hosts a low-luminosity AGN in its nucleus. Additionally, 
there is abundant evidence pointing to vigorous star formation tak- 
ing place in the inner ~ 20" x 30" region of this galaxy (e.g., 
;|Poggelll989l ; iTelesco et all 1 19931 ; iDavies et alj 
2005). This is also a typical 'single-peak' 



Puxley etalj|l988 



19971 ; IJogee et all 



type 

galaxy, with the CO molecular gas being concentrated in a nuclear 
disc of ~ 10" radius and an unresolved co mpact core at the nucleus 
(e.g JSofue et alj|2003l ; IJogee et al.120051) . Previous low-resolution 
(1.6") images of H 2 2.1 2 /im emission in the inner 10" of this 
galaxy were analysed bv lDavies et alj Jl997l) . They found a strong 
NIR continuum peak in the nucleus, but the H 2 emission has peaks 
distributed in what seems to be an edge-on ring structure around the 
nucleus. Their analysis suggests that while active star formation is 
taking place in the circumnuclear ring, the star form ation in its nu- 
cleus has come to an en d (see also ljogee et alj20"05l) . Very recently, 
iRosenberg et al.l ( f2012h presented non-AO H-, J- and AT-band SIN- 
FONI data of the inner 8" x 8" of NGC 4536. These authors con- 
firmed the presence of a circumnuclear ring in NGC 4536, traced 
mainly by Bry and [Fen] 1.26 jim emission, while the H 2 2.12 /jm 
distribution is more centrally peaked. 

High- and low-resolution data of the nuclear region of 
NGC 4536 were obtained with SINFONI, covering the inner 3"x3" 
and 8" x 8", respectively. Intense Bry and H 2 emission lines 
were detected in the spectra of this galaxy, together with weaker 
He 1 2.06 yum emission. The flux distribution maps of H 2 2.12 jim, 
Bry and Hei derived from the SINFONI data are shown in Fig. [6] 
The left panels correspond to the maps derived from the high- 
resolution data (HR maps) and the right panels show the ones de- 
rived from the low-resolution data (LR maps). 

It can be seen in Fig. [6] that the circumnuclear molecular and 
atomic gas of NGC 4536 exhibit very different morphologies. The 
Bry emission is mainly localised in a circumnuclear ring ~ 7.5" 
(540 pc) in diameter, oriented along the photometric semi-major 
axis of the galaxy. The ratio between the semi-major and minor 
axis of the ring (~ 2.8) is consistent with a circular ring with the 
same inclination as the galaxy (i = 69°, Table [2} projected onto 
the sky. The Bry emission is not uniformly distributed along the 
ring but shows a knotty appearance, with the strongest emission 
located towards the NW. In contrast, the H 2 emission is more cen- 
trally concentrated, with the bulk of the emission restricted to the 
inner ~ 2" of the galaxy. The detailed structure of this region can 
be seen in the HR H 2 distribution map, which shows that the H 2 
emission does not peak at the position of the active nucleus but 
seems to encircle it. The nuclear region is connected to the circum- 
nuclear ring by two ridges of weak H 2 emission extending towards 
the south and the north, giving a global S-shaped appearance to the 
H 2 distribution. In general, the morphologies of the H 2 and Bry 
emission-line s showed by the LR m aps are consistent with those 
presented by IRosenberg et alj J2012h . The presence of strong H 2 
emission in the nucleus of NGC 4536 and the absence of strong 
hydroge n recombination line s (i.e., Bry) support the scenario pro- 
posed bv lDavies et alj dl997l) . in which no current star formation is 
taking place in the nuclear region of this galaxy (see Section |4~2l for 
a detailed discussion). 

NGC 4536 is the only galaxy in the sample for which 
He 1 2.06 yum line emission was observed. The morphology of the 
He 1 line-emitting gas (lower-right panel of Fig. |6]l has an almost 
perfect one-to-one correspondence with the Bry emission, indicat- 
ing that both Bry and He 1 are tracing the same structures. 

While previous low-resolution observations sug gested the 
presence of a circumnuclear ring in this galaxy (e.g., IVila et all 
ll990l ; lDavies et alJl997l ; lLaine et alj2006h . the high spatial resolu- 
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Figure 7. Flux distribution map of the H2 2.12 fim emission line observed 
in NGC 4569. The dashed line indicates the position of the major axis of 
the galaxy and the dotted line corresponds to the stellar bar position angle. 
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Figure 8. Flux distribution map of the H2 2.12 fim emission line observed 
in NGC 4579. The dashed line indicates the position of the major axis of 
the galaxy and the dotted line corresponds to the stellar bar position angle. 



tion provided by SINFONI gives a clear view of the circumnuclear 
ring present in the inner ~ 8" of NGC 4536. A more detailed study 
of the ionized gas ring of this galaxy is presented in Section |4.2.1| 

3.5 NGC 4569 

NGC 4569 (M90) is a member of the Virgo cluster classified 
as a SAB(rs)ab galaxy. Its nucleus was classified as a 'transi- 
tion object' or weak-[Oi] LINER by iHo et all i ll 9971) . although 
the mechanism producing the nuclear activity is still unclear, with 
several lines of evidence indicat ing a strong contribution from 



199S 




2001 





Ho et al J 1 19971: iMaoz et al 



..Keel||199L.. . 

Barth & ShieldlbOOd; lAlonso-Herrero et alJl2000l : lHo et al 
Gabel & BruhweilefeOO^ This galaxy has a bright, point 



like nucleus at optical and UV bands that dominates the emis- 
sion. The HST V - I map show that while the circumnuclear 
regio n is dusty, the nu cleus itself is apparently unobscured by 
dust jPogge et alj|2000l) . NGC 4569 is considered an 'anemic spi- 
ral' because of its Hi deficiency (it has ~ 1/10 of the atomic 
gas that similar field g alaxies have) and small extent compared 



to its optical size (e.g.. | van den Berghlll976l : ICavatte et alj| 19941 ; 



iKoopmann & Kenney||2004[). Its molecular CO g as has been ex 
tensively studied (e. g., Kermev & Yound [l986l; I Sakamoto et al] 



1999b a; Sofue et al. 2003 : iHelfer et al J 12003; Uogee et al J 12005 
Nakanishi et al.l2005l : lBoone et alj2007li201 llJwilson et al.l2009T) . 

The CO emission is highly concentrated in a compact bar-like re- 
gion (of ~ 17" x 6" extending along P. A. as 24°) with three off- 
centre peaks, one located at ~ 1" towards the north-e ast of the nu- 
cleus and two symmetrical peaks separated by ~ 10" l lBoone et alj 
120071) . 

The SINFONI observations of NGC 4569 were done in 
the high-resolution configuration, covering a total FOV of 3" x 
3". The A"-band emission line spectrum of this galaxy is domi- 
nated by H 2 transitions, with no signs of emission from ionized 
gas. This is consistent with previous long-slit observations where 
no Pa/? and/or Bry emission was detected in the nucleus (e.g., 
lAlonso-Herrero et al .l2000l : lRhee & Larkinl2005h . Fig.|7]shows the 



H2 flux-distribution map derived from H2 2.12 fan. It can be seen 
that the H 2 emission line gas distribution is very irregular, with 
most of the emission located off-nucleus on the eastern part of the 
FOV. There is a rapid drop in the signal-to-noise as one goes from 
the eastern-southern regions towards the north, where the intensity 
of the H 2 line falls below the continuum noise level. Although the 
SINFONI data presented here is mapping only a small part of the 
molecular bar observed in NGC 4569, it is possible to say that the 
general H 2 distribution displayed by this gala xy is consistent wit h 
the CO(1-0) and CO(2-l) maps reported by iBoone et"al] d2007t) . 
in the sense that the stronger emission is located towards the east 
and south. However, the H 2 emission peaks closer to the nucleus 
than the CO emission, suggesting a strong gradient in the molec- 
ular gas temperature going from the centre to the outer regions of 
the galaxy. 

3.6 NGC 4579 

NGC 4579 is an SAB (rs)b galaxy with a LINER/Seyfert 1.9 nu- 
cleus jHo et aljfl997l) . As in the case of NGC 4569, this galaxy 
is considered an 'anemic spiral' of the Virgo cluster. HST UV, 
[Om] and Ha images show a bright nuclear point source sur- 
rounded by filamentary structures, with the brightest emission 
located at distanc es < 2" towards the NE dMaoz etal.lll995l : 
IPogge et alj|2000h. High-reso l ution CO intensity maps presented 

showed that most of the cold 



Garcia-Burillo et al 



gh-resolution 1 
T d2005Ll2009h 



by 

molecular gas in NGC 4579 is located in a two spiral arm structure 
that can be traced from r ~ 10" down to a distance of 2" from the 
nucleus. The spiral arms lie at the leading ed ges of the ~ 12 kpc- 
diameter bar identified in K-bsmd images (e.g jKnapen et al.l2003l) 
and a re well correlated w ith dust lanes revealed by V - I HST im- 
ages dPoggeetal.ll2000h . especially the arm located towards the 
north. Some emission from CO is observed at smaller distances 
from the centre but not on the position of the AGN it self. 

The SINFONI Zf-band spectra of the inner 3" x 3" of 
NGC 4579 display strong H 2 emission-lines and no evidence of 
ionized gas. Two main structures are seen in the H 2 emission-line 



© 0000 RAS, MNRAS 000, 000-000 



High resolution 
Flux [1e-18erg s 1 cm" 2 ] 



Molecular gas in the centre of galaxies - I 9 
NGC 4536 

Low resolution 
Flux [1e-18erg s 1 cm" 2 ] 



o.o 



2.3 



4.7 



7.0 



o 

CD 

o U 

1 

CO 



* 



-1 

arcsec 



Flux [1e-18erg s 1 cm" 2 ] 



0.0 



3.5 



7.0 10.5 



o 

CD 

0) n 

o U 

i 

CO 




-1 
arcsec 





0.0 


11.3 


22.5 


33.8 




-4-2 2 
arcsec 

Flux [1 e-1 7erg s 1 cm" 2 ] 



0.0 



2.4 



4.8 



7.2 



CD „ 

cn n 
o u 

i 

CO 



+ 9 



V 



-4-2 2 
arcsec 

Flux [1e-18erg s 1 cm" 2 ] 



0.0 



13.6 27.3 



40.9 



o 

CD 
C/3 
O 
i— 
CO 



m 

Hel 



-4-2 2 
arcsec 

Figure 6. Flux distribution maps of the Ht 2.12 jum (upper panels), Bry (middle panels) and He 1 2.06 fim (lower panel) emission lines observed in NGC 4536. 
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map shown in Fig. [8] This emission is dominated by a symmetric 
nuclear component, slightly displaced towards the S with respect to 
the maximum of the continuum. Towards the NE, at a distance of 
~ 1" from the nucleus, a weaker component is present. The overall 
shape of the H 2 emission resembles the one revealed by Ha and 
[Om] narrow-band images, although the latter show much more 
comp lex structures resol ved into individual knots (see figure Id and 
3a of IPogge et al"1l200d) . Additionally, there is no correspondence 
between the H 2 emission presented here with dust or CO emission, 
which is mostly found at larger distances from the nucleus. 



4 PHYSICAL CONDITIONS OF THE NUCLEAR 
REGIONS 

Rotational and vibrational H 2 transitions are important cooling 
channels for warm molecular gas, observed as a series of IR emis- 
sion lines. H 2 emission lines are commonly observed in the NIR 
spectra of galaxies harbouring energetic processes related to mas- 
sive star formation, shocks or AGNs, and have proved to be impor- 
tant diagnostic tools in all molecular media. 

Most of the galaxies in our sample show signs of AGN activ- 
ity in their centres and, at the same time, regions where star for- 
mation is taking place. Therefore, the data presented here give us 
an excellent opportunity to study the physical conditions of the H 2 
emission-line gas in these two distinct environments. In order to 
achieve this we extracted spatially integrated ID spectra from the 
nuclear and circumnuclear regions of the galaxies. The size and lo- 
cation of these regions vary from galaxy to galaxy depending on 
their particular morphology and are indicated in the first column of 
Table [3] Many emission lines were detected in the SINFONI spec- 
tra of the sample galaxies, including Bry and H 2 1.96, 2.03, 2.12, 
2.22 and 2.25 /jm. As described below, these lines provided us with 
a variety of tracers of different physical quantities. The emission- 
line fluxes of the individual lines were obtained by fitting a single 
Gaussian comp onent plus a linear continuu m to the spectra. The 
'Liner' routine JPogge & De Robertid[l99~3r) was used for this pur- 
pose. The measured emission-line fluxes can be found in Table [3] 
The precise position and size of the apertures used to extract the ID 
spectra is given in Table lAll of Appendix lAl 

4.1 Excitation mechanisms and temperature of the H 2 
emission-line gas 

Several mechanisms have been proposed for the excitation of the 
ground state ro- vibrational levels of H 2 . These can be separated 
into two main groups: thermal and non-thermal processes. In ther- 
mal processes, the H 2 molecule is excited via inelastic collisions 
with atoms or molecules in warm (T~ 500-2000 K) gas, heated ei- 
ther by shocks or UV/X-ray radiation in high density r egions (e.g., 
Lepp & McCravl ll983l:JSterriberg & Dalgarnd Il989l : iBrand et alj 
19891 : iDraine & Woodslll99Cl : iMalonev et al.lll996^ In this case, 
to a first approximation, the H 2 line-emitting gas is in local ther- 
mal equilibrium (LTE) and the spectrum is therefore characterized 
by a single excitation temperature. Since H 2 is rapidly dissociated 
by thermal collisions above ~ 4000 K, the excitation tempera- 
ture of the thermal component must be lower than this. In non- 
thermal processes, the NIR ro-vibrational lines are observed due 
to a radiative cascade through the excited vibrational and rota- 
tional levels of the ground electronic state. The excitation of the 
electronic state can be produced either by near-UV (11-13.6 eV) 
pumping to the electronically excited Lyman and Werner bands 



(e.g.. lBlack& van Dishoeckll987l) or by a collision with a fast elec- 
tron produced by X-ray i onization to singlet electronic states (e.g., 
iGredel & Dalgarnolll995l) . 

Each of the excitation mechanisms mentioned above produces 
different H 2 spectr a, which are often observed together within a 
single object (e.g 



2005 



Moui'i 



19941: lyeilleux et alJI 19971: iDavies et alj 



1997L 120051: lOuillen et al l 1 19991: iRodrfguez-Ardila et alj 



Riffel et alj 



2006 



2008, l2010l : IZuther et al 



2004 



2007 



Storchi-Bergmann et alj 20091) . The relative intensity of the 



H 2 emission lines can be used to infer the dominant pro- 
cess responsible for the H 2 emission. In particular, the H 2 
2-lS(l) 2.25 /jm/l-0S(l) 2.12 yim line ratio is an excellent 
discriminator between non-thermal and thermal processes, while 

1- 0S(2) 2.03 ^m/l-0S(0) 2.22 yum can be used to differentiate 
between the rmal UV excitation and X-ray/shock excitation 
jMouri|[l994l) . Additionally, these line ratios are sensitive to the 
vibrational temperature (T v ;t,) and rotation al temperature (T rot ) 
of the emission-line gas, respectively (e.g.. IScoville et all 1 19821 : 
iReunanen et alj 120021) . The excitation temperature (T vib Or 1 ro t , 
depending on the transitions involved) will depend on the relative 
importance of the radiative and collisional processes. In the case of 
gas in LTE the ro-vibrational level populations are characterized by 
a thermalised Boltzman population and the excitation temperature 
corresponds to the kinetic temperature of the gas . Assu ming 
the computed A-Einstein coefficients of iTurner et alj l ll977l) . the 
relation between the excitation temperature of the gas and the H 2 

2- lS(l)/l-0S(l) and 1-0S(2)/1-0S(0) emission line ratios can be 
expressed as: 



5600 



In (1.355 x^£fii) 



-1113 



ln (0 .323 x^i)' 

'1-05(0) 



(1) 



(2) 



Fig.|9]shows the H 2 2-lS(l)/l-0S(l) vs. 1-0S(2)/1-0S(0) line-ratio 
diagnostic diagram for the selected regions of the galaxies in our 
sample, together with diverse key excitation models taken from the 
literature. Each filled symbol represents one particular galaxy and 
the associated numbers indicate the different regions specified in 
the first column of TablefS] The right and upper axes show the val- 
ues of T vib and T rot , respectively, derived from Eqs. (TJ and {2jl as a 
function of the two line ratios. Table[3]lists the temperature values 
obtained for the different regions of the galaxies. From this figure 
it is clear that there is no single excitation mechanism acting in the 
nuclear and circumnuclear regions of the galaxies, regardless of 
their type of activity. In general, each region is characterized by a 
T rot lower than its T vib , consistent with some contribution from non- 
thermal processes. The values of temperature derived vary from 
~ 700 up to 4000 K, with a tendency for the objects to be located 
near the pure thermal excitation curve for the range of tempera- 
tures ~ 1000 - 3000 K. These high temperatures correspond to the 
regio n on the diagram w here X-ray/shock thermal excitation domi- 
nates. |Dorreral]|20l2) const ructed a grid of phot oionization mod- 
els using the code cloudy/08 dFerland et alJll998o . which includes 
both collisional an d radiative interactions of H 2 and its environment 
dShaw et al]|2005l) . The input parameters are a spectral energy dis- 
tribution of the ionizing source similar to that observed in typical 
AGNs, three values of metallicity (0.5, 1 and 2 solar metallicy), 
an electronic density n e = 10 4 cm" 3 and an ionization parameter 
U in the -4.0 < logC < -1 range. The resulting H 2 line ratios 
span a large range of values, H 2 2.25/2.12 /jm ~ 0.05 - 0.5 and 
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2.03/2.22 /im ~ 1.1 - 2.0, reproducing most the observed H2 line 
ratios of the galaxies in our sample. 

While most of the galaxies in our sample are characterized by 
a value of H2 2.03/2.22 yum higher than 1, this ratio is lower in 
the case of NGC 4536. The H 2 line ratios of this galaxy occupy a 
region on the diagram dominated by a mix of UV thermal and non- 
thermal models, suggesting that UV radiation plays a key role in the 
excitation of the H 2 in this galaxy, with no significant contribution 
from shocks and/or X-ray radiation. The combination of low T rol 
and high T vib make this galaxy unique among the objects presented 
here, placing it in a region of the diagram typically occupied by 
photodissociation regions. 

In summary, the analysis presented in this section suggest that 
thermal processes are primarily responsible for the excitation of the 
H2 gas in the nuclear and circumnuclear regions of the galaxies in 
our sample, with some minor contribution from non-thermal pro- 
cesses (i.e., UV fluorescence). 



4.2 Ionized gas and the H 2 /Bry emission line ratio 

The Bry recombination line was observed in half of the galaxies 
of the sample and, even though the number of objects is low, we 
can appreciate very different morphologies displayed by this line. 
In two cases (NGC 3351 and NGC 4536) the flux distribution maps 
of Bry show a clear ring-like morphology, while in the third case 
(NGC 3627) the bulk of emission is restricted to a relatively com- 
pact nuclear region. 

In general, the Bry distribution of the galaxies in the sample 
does not follow their H 2 emission. A clear deficit of Bry emission 



relative to the H 2 emission is observed in the nuclei of NGC 3351 
and NGC 4536. In order to quantify this, we included in Table [3] 
the values of the H 2 2.12 jimfBiy emission line ratio measured for 
the different regions of the galaxies. These values vary between 
0.45 and 5.9, in agreement wit h the usual values found for star- 
burst galaxies and AGNs (e.g., Moorwood & Olival 1 19881 1990l. 



Larkin et al ll998l ; lRodrfguez-Ardila et alj200 4. 2005; Riffe l et alj 
20101) . The expected line ratio for starburst-dominated galaxies is 
typically low (< 0.6), while LINER nuclei show the highest ra- 
tios, usually > 2. Composite and Seyfert nuclei are usually char- 
acterized by a H 2 /Bry ratio of order unity. The lowest of our val- 
ues are the ones measured in the circumnuclear regions and agree 
with the expected values for pure starburst emission. On the other 
hand, the highest values are the ones measured in the nucleus of 
the galaxies, which are in the range usually occupied by LINER 
nuclei. While this is consistent with the nuclear classification of 
NGC 3627 and NGC 4536, the galaxy NGC 3351 is known to be 
a non-AGN. However, the high value of H 2 /Bry displayed by the 
nucleus of this galaxy can be explained as a consequence of the rel- 
ative old stellar population that characterizes its innermost regions 
(see below), in which the presence young OB stars responsible for 
the Bry emission has declined and the number of SN explosions, 
which heat the gas via shocks and produce the H 2 emission, has 
increased. Although a detailed interpretation of the variation of the 
H 2 /Bry from pure-starburst emission to Seyfert and LINER nuclei 
is non-trivial, one appealing scenario is that the different line ratios 
reflect different excitation mechanisms, in the sense that they de- 
scribe a transition from purely ionizing radiation powered by star 
formation to pure shock excitation driven by SNRs. The H 2 /Bry 
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Table 3. Emission line fluxes (in units of 10 16 erg cm 2 s ), temperatures, emission-line ratios and molecular masses measured for different regions of the 
galaxies of the sample. The errors correspond to 3tr. 



Galaxy 


H 2 1-0S(2) 


H 2 1-0S(1) 


H 2 1-0S(0) 


H 2 2-lS(l) 


Bry 


Tvib 


T rM 


H 2 2 .1 2 fim 


M 


M co ld 


2.03 yum 


2.12 fim 


2.22 fim 


2.25 fim 


2.16 fjm 


[K] 


[K] 


Bry 


[Mo] 


[10 6 M o ] 


NGC 3351 






















1 - Nucleus 


3.9 ± 1.3 


7.5 ± 1.6 


3.4 ± 1.2 


1.2 ±0.9 


1.3 ±0.9 


2600 


1125 


5.9 


3.9 


2.8 


2 - W cloud 


0.7 ±0.1 


1.2 ±0.2 


0.5 ±0.2 


0.2 ±0.1 


- 


2800 


1500 


- 


0.6 


0.4 


3 - 3" X 3" 


21.4 ±5.8 


16.7 ±6.7 


11.9 ±6.7 


< 6.0 


- 


<4300 


2000 


- 


8.7 


6.2 


4 - 8" x 8" 


117 ±55 


151 ±26 


54.3 ± 28.7 


37.5 ±28.1 


150 ± 25 


3300 


3000 


1.0 


79.0 


57.0 


5 - Ring 


- 


89.0 ± 22.4 


- 


- 


160 ± 28 


- 


- 


0.56 


46.5 


33.3 


NGC 3627 






















1 - Nucleus 


18.2 ±6.9 


43.1 ±6.0 


15.0 ±5.7 


6.2 ±5.7 


20.6 ± 6.6 


2500 


1200 


2.1 


22.3 


16.0 


2 - 3" X 3" 


41.8 ±22.9 


108 ± 18 


32.1 ± 17.4 


< 17.0 


- 


< 2600 


1300 


- 


55.6 


40.0 


3 - N cloud 


4.7 ± 0.9 


11.7 ±0.5 


3.0 ±0.6 


1.3 ±0.6 


- 


2207 


1680 


- 


6.0 


4.4 


4 - S cloud 


6.0 ± 0.8 


15.2 ±0.9 


4.2 ±0.9 


1.4 ± 1.0 


- 


2091 


1388 


- 


7.S 


5.6 


5-8" X 15" 


- 


268 ± 62 


- 


- 


- 


- 


- 


- 


138 


100 


NGC 4501 






















1 - Nucleus 


- 


7.6 ±3.9 


5.8 ±4.5 


- 


- 


- 


- 


- 


10.5 


7.5 


2 - Ring 


48.1 ±27.0 


50.1 ± 12.9 


27.1 ± 18.0 


< 18.0 


- 


< 4200 


2000 


- 


69.4 


49.9 


3 - 3" X 3" 


- 


65.4 ± 20.3 


- 


- 


- 


- 


- 


- 


90.6 


65.0 


NGC 4536 






















1 - Nucleus 


11.5 ±6.2 


32.6 ±5.6 


17.1 ±5.9 


9.0 ±5.9 


17.0 ±5.2 


3520 


730 


1.9 


37.0 


26.6 


2 - Ring 


92.2 ± 26.0 


240 ± 29 


105 ± 22 


88 ±25 


532 ± 29 


4300 


885 


0.45 


272.6 


196.0 


3 - 8" X 8" 


140 ± 59 


331 ±34 


151 ±36 


114 ±36 


567 ± 33 


4000 


920 


0.6 


375.9 


270.0 


NGC 4569 






















1 - 3" x 3" 


52.2 ±21.2 


95.4 ± 18.8 


42.5 ± 19.5 


24.8 ± 19.5 




3400 


1200 




74.1 


53.3 


NGC 4579 






















1 - Nucleus 


29.5 ±4.5 


46.8 ±7.5 














64.8 


46.6 


2 - NE cloud 


10.2 ±2.3 


19.0 ±2.4 


5.3 ±2.5 


3.0 ±2.4 




2600 


2300 




26.4 


19.0 


3 - 3" X 3" 


61.8 ±25.8 


82.9 ±21.3 


34.9 ± 19.7 


< 22.5 




3500 


> 2000 




114.9 


82.6 



derived for the galaxies together with the analysis of the excitation 
mechanisms presented in the previous section are compatible with 
this picture: the lowest H 2 /Bry ratios are displayed by the outer 
regions of NGC 4536, which are the ones with the strongest non- 
thermal contribution (UV fluorescence) to the excitation of the H 2 
molecule, while the highest H 2 /Bry ratios are measured for regions 
with a strong influence from thermal processes (e.g., shock heat- 
ing). 



4.2.1 Star-forming rings in NGC 3351 and NGC 4536 

The SINFONI data presented here allowed us to spatially resolve 
the innermost regions of NGC 3351 and NGC 4536, both of which 
show circumnuclear rings traced by ionized gas (i.e. Bry and/or 
He i emission lines). This type of ring is frequently found in early 
to intermediate Hubble-type barred galaxies, and is usually associ- 
ated with inflowing gas accumulating near the inner Lindblad reso- 
nances of bars (see lButa & CombesHl996l and references therein). 
Photoionization by OB stars is the most likely source of the neb- 
ular emission, although shocks produced by stellar winds and su- 
pernovae (related to the same stellar population) may play an ad- 
ditional role. In addition, high energy photons from AGN can 
also ionize the ISM and contribute to the Bry emission. Although 
NGC 4536 is known to harbour a low-luminosity AGN, the spatial 
distribution of its Bry emission (Fig.[6j indicates that this contribu- 
tion is relatively low and restricted to the innermost regions, with 
no or negligible influence in the circumnuclear ring observed in 
this galaxy. Therefore, outside the nucleus we consider Bry as a 
primary tracer of young stellar populations. 

It is possible to estimate upper limits of the age of the star 



formation occurring in these galaxies by measuring the equivalent 
width of the Bry line, Wb, t . For this, we extracted two spectra for 
each galaxy, one integrating the emission of the inner ~ 1" of the 
galaxies and a second one encompassing the emission coming from 
the ring. Note that in the case of NGC 3351 only part of the circum- 
nuclear ring is covered by the field of view of our data (Section [3~TT l. 
We measured W Bry of 0.3 and 2.1 A for NGC 3351 and 0.8 and 
8.3 A for NGC 4536, with the higher values corresponding to the 
ring spectra. Using the Starburst99 models jLeitherer et al. I [19991) 
for an instantaneous burst, different slopes for the initial mass func- 
tion (2.35 and 3.33), upper mass limits of 30 and 100 Mq and solar 
metallicity, the measured values of W B , T predict ages of the stellar 
populations of 10.5 and 6.9 Myr for NGC 335 1 and 7.2 and 6.5 Myr 
for NGC 4536, with the younger populations located in the cir- 
cumnuclear rings. As mentioned above, the age values correspond 
to upper limits, since we have not taken into account the possible 
contribution from the AGN and/or underlying (old) stellar popu- 
lation to the continuum. This additional component increases the 
continuum level, producing lower values of Wsry than the true ones 
and, hence, larger ages than the expected for pure star-forming re- 
gions. An estimation of the contribution of the AGN to the con- 
tinuum can be made based on the values derived for the equiva- 
lent width of the CO(2-0) absorp tion line at 2.29 fim, Wcop-O), as 
proposed by iDavies et alj d2007l) . These authors showed that the 
intrinsic Wcoa-0) of any stellar population is typically ~ 12 A 
(with an uncertainty of ±20 per cent), independent of the star for- 
mation history and age. Therefore, for a stellar continuum diluted 
by additional non-stellar emission, the fraction of the stellar light 
can be estimated by the ratio of the observed and intrinsic equiv- 
alent widths of the CO feature. In order to check for a possible 
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AGN contribution to the continuum in NGC 4536 we measured 
the W C o(2-0) m tne nuclear and ring spectra of this galaxy. The 
values are W C o(2-0) = (13.6 ± 0.8) A for the nuclear region and 
(1 4.0 ± 0.7) A for the ring. These values are in the range derived 
by iDavies et alj ( |2007| ). suggesting that the contribution from the 
AGN to the observed continuum in NGC 4536 is negligible. 

Alternatively, we could assume a scenario where continuous 
star formation is taking place. In this case, the Starburst99 models 
predict ages of more than 60 Myr for the measured W Bry . Such 
old ages are improbable, since the star formation would need to be 
sustained over long periods of time in very localized regions, in 
which turbulence and heating from SNe will disrupt the molecular 
clouds after ~ 2 - 30 Myr, inhibiting further star formation (e.g., 
iBlitz et alj|2007l) . 

In addition to Bry, the spectra of NGC 4536 also display 
He i 2.06 yum emission along the circumnuclear ring (Fig. [6^. The 
ionization energy required to produce Hei recombination lines is 
24.6 eV, almost a factor of two higher than that of the hydrogen 
(13.6 eV). Therefore, Hei emission lines are expected to arise in 
the vicinity of very hot massive stars (most massive OB stars) and 
AGNs, where enough high-energy photons are available to ionize 
the He i gas. Since the hottest stars will vanish fastest, the He i/Bry 
emission-line ratio can give us a qualitative estimation of the rela- 
tive ages of the young star clusters located along the circumnuclear 
ring observed in NGC 4536. We measured a value of He i/Bry ^ 0.5 
along the entire ring, showing no significant variations and/or gra- 
dients. Hence, the young star clusters along the ring of this galaxy 
have the same ages, suggesting a scenario where the star formation 
started at the same time over the entire structure. It must be kept in 
mind that this is only a qualitative interpretation, since a detailed 
assessment of the He i/Bry ratio is complicated by additional de- 
pendences such as the geometry of th e ionized gas, its density, dust 
conte nt and helium abundances (e.g., IShieldsll993l ; lLumsden etafl 
1200 lh . 



5 NUCLEAR MASS INVENTORY 
5.1 Warm molecular gas mass 

The total mass of the warm molecular hydrogen, M warm , can be de- 
rived from the observed flux of the H 2 1-0S(1) 2.12 yum emission 
line. Assuming a temperature T = 2000 K, the H2 1-0S(1) transi- 
tion probability As(i) = 3.47 x 10~ 7 s~\ and a population fraction 
in the (y, J) = (1,3) level f {l 3) = 0.0122, the relation between the 
warm molecular mass and the H 2 1-0S(1) emission line flux is given 



by: 



13 f— ) 2 (" 



Fi_ 



0S( I ) 



10" 



(3) 



where D is the distance to the g alaxy and A 2 2 is the extinction at 
2.2 fj.m (e.g., IScoville et alj 19821) . Unfortunately, the data presented 
here do not allow us to determine the extinction affecting the inner 
regions of the galaxies in our sample and therefore it is only pos- 
sible to determine a lower limit for M waral , assuming a value of 
A 2 2 = 0. Taking into account that the absorption at the A"-band is 
only » 10 per cent of that at visible wavele ngths dCardelli et al] 
1 19891) . and a typical extinction A v ~ 1.2 mag dHo et alJI 19971) . cor- 
recting for internal extinction will not increase the derived masses 
by more than 12 per cent (or a factor of 2 for an extreme case of 
A v ~ 7.7). Table[3]lists the integrated masses of warm molecular 
gas derived for the different regions of the sample galaxies studied. 



These cover a wide range of values, ranging from ~ 9 M© in the 
inner 75 pc of NGC 3351 up to ~ 300 M Q in the inner 280 pc of 
NGC 4536. 



5.2 Cold molecular gas mass 

In principle, the H 2 1-0S(1) 2.12 /jm emission does not necessarily 
reflect the cold gas distribution, since it also depends on an external 
source of energy capable of exciting the H 2 ro-vibrational levels, 
such as UV photons from young OB stars, shocks driven by SNRs 
or by outflows or winds from the active nucleus (see Sect. 14.1b . 
Historically, the favoured proxy for estimating the amount of cold 
H 2 gas has been CO emission. However, it has been suggested 
that an estimate of the total content of cold H 2 gas can be ob- 
tained from the H 2 2.12 yum line luminosity (Li_osu))- This method 
has the advantage of being able to probe the inner-most regions 
of near by objects with higher spatia l resolution than CO obser- 
vations. iMuller Sanchez et alj j2006t) derived a mean conversion 
factor of 4000 Mq/ L© (with a standard deviation of a factor of 
~ 2) from the relation between the H 2 2.12 luminosity and 
the H 2 mass traced by CO luminosity in a sample of sixteen lu- 
minous and ultraluminous infrared galaxiefl Taking into account 
Eq. 10, these values yield cold-to-warm H 2 mass ratios in the range 
of M cold /M warm =a (1 -5)x 10 6 . On the other hand. lDale et alj J2005h 
studied a large sample of active and star forming galaxies and found 
cold-to-warm H 2 mass ratios covering a much larger range of val- 
ues, M coM /M wal . m a 10". 

In order to investigate the relation between the H 2 2.12 yum 
emission and the cold H 2 mass, we compiled from the literature 
values of M coU derived from CO observations and H 2 2.12 yum lu- 
minosities for a large number of galaxies, covering a wide range 
of luminosities, morphological types and nuclear activity. Panel (a) 
of Fig.fTOlshows the values of M co i d derived from CO observations 
as a function of the integrated area (i.e., size of the aperture used 
to integrate the CO flux). The relation between the cold H 2 gas 
mass and the integrated area, A, can be described as a power-law 
M coM = a A 7 , with the best-fitting parameters a = (1.0 ± 0.3) x 10 6 
and y = 0.39 ± 0.06. Many factors contribute to the large dispersion 
observed in this relation (cr(log(M)) = 0.5), including uncertainties 
in the CO-to-H 2 conversion factor (which likely differs from galaxy 
to galaxy), the inhomogeneity in the type of galaxies included, and 
the intrinsic irregularity in the gas distribution of the galaxies (e.g. 
presence of bars, rings). Similarly, panel (b) of Fig. 1 101 shows the 
distribution, as a function of the integrated area, of the H 2 2.12 yum 
luminosities compiled from the literature. These follow a similar 
trend as the cold H 2 masses, with a value of y = 0.49 ± 0.07 con- 
sistent at lcr with the one derived for the Mc id distribution. This 



One outlier, NGC 6240, with a corresponding factor of only 400 Mq/ Lq 

was taken out of the sample. 

3 Reference s for the C O data: ISakamo to et al 



IWalter et alj feOOll |200 j): ISchinnerer et alj j2002E 



ll999bfl 



Garcfa-Burillo et al 



|2003|) ; lGao & Solomonlj2004l);ICombes et alJfe004l2009l) : IOnodera et al 

Ijogeeet alj 12~ ~ 



20041) : iGara a-Barreto et al. 1 20051): Uogee et alj 120051): Boone et al 
llono et alj 120071): Ifficks et all J2009h: llsraell J2009I): Izhu et al 
Casasola et alj J2010l |201 lh T lMartm et alj j2O10l) : iMeier et al 




Olsson et al 



Salome et al 



Dale et al 



20071 l2Q10l) ~ lAlatalo et alj 1201 ll) : ICrocker et al 
I20T l|). References for H 2 2.12 11m line lumf 



20041) 



Rodriguez-Ardila et al 



Miiller Sanchez et al 
iNowak et alj 120081) 



(2006 



Riffel et al 



IStorcM-Bergmann et alj J2009): Frie drich et al 



|Daviesetal 
Zuther et alj <2007l): iBbker et al 
l2OO8tl2O10l):lHicks et alj 120091) 



2010) 
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10 2 10 3 10 4 10 5 10 6 10 7 10 B 10 9 

area/ pc a 

Figure 10. Panel (a): Cold H2 gas masses derived from CO observations as 
a function of the integrated area. The solid line corresponds to the power- 
law that best reproduces the data. Panel (b): H2 2.12 /im luminosities as a 
function of the integrated area. The solid line corresponds to a power-law 
fit to these data. Panel (c): Cold H2 gas masses derived from H2 2.12 /im 
luminosities using a rescaling factor ji = 1174 Mq/ Lq. The solid line 
corresponds to the same power-law as in panel (a), while the dashed line 
corresponds to the relation that would be obtained using a factor of /? = 
4000 Mq/ Lq (see the text for details). Red squares correspond to values of 
the galaxies in our sample. 



is in agreement with the previous suggestion of a proportionality 
between the warm and cold molecular gas masses and gives us an 
opportunity to obtain a rough estimation of the factor relating these 
two quantities. Assuming that the two relations have the same slope 
(y = 0.39), the proportionality factor can be estimated by minimiz- 
ing the residuals of the rescaled luminosities,/? x L|_os(u, with re- 
spect to the power-law that best describes the distribution of cold 
H 2 masses derived from the CO measurements. Panel (c) of Fig.fTOl 
shows the values of Lh 2 compiled from the literature rescaled by the 
factor {3 = 1 174 Mq/ Lq obtained in this way. We also included the 
re sult obtained using the r escali ng factor of 4000 Mq/ Lq proposed 
bv lMuller Sanchez et all | |2006|) (dashed line). Note that, if we had 
used that factor, we would had overestimated the masses as com- 



pared to the ones derived from the CO observations. Therefore, the 
factor of 1174 seems to give more representative cold H 2 masses 
for a larger range of galaxy types. This factor and its corresponding 
error (cr\og(J3) = 0.35) can be translated to a range of M co i,i/M walm 
ratio of ^ (0.3 - 1.6) x 1 6 , which is consistent with the one de- 
rived bv lDale et al .1 1|2005[) and margina lly consistent with the one 
inferred by Miiller Sanchez et all d2006h . 

From the analysis described above, we can conclude that a 
rough estimate of the cold H 2 gas mass can be obtained from the 
integrated H2 2.12 //m luminosity using the relation 



Me, 



M„ 



« 1174x 



Li _i 



0S( 1 ) 



(4) 



Taking this relation into account, we calculated the total (cold) 
H 2 gas content for the galaxies in our sample. The observed 
H2 2.12 /mi line fluxes and corresponding cold LL masses of the 
nuclear and circumnuclear regions of the sample galaxies are listed 
in Table [3] We found masses between 4 x 10 5 and 2.7 x 10 s Mq 
within regions with radii of up to 280 pc. In order to place these 
results within a larger context, we included in the lower panel of 
Fig.[l0]the cold H2 masses derived from the H2 2.12 jim luminosi- 
ties measured for the largest regions in each galaxy. The masses 
obtained for these galaxies agree very well with the 'global' trend 
followed by the large pool of masses taken from the literature. The 
only galaxy that seems to deviate from this trend is NGC 3351, 
for which the inner regions show a lower cold-to-warm mass ratio 
than the rest of the sample. However, this ratio goes up when it is 
calculated within a larger radius region (i.e., 8" x 8")- 

The values of M co id derived from the LL luminosities can 
be compared with values derived from CO observations. For 
this, we also included in the upper panel of Fig. \W\ val- 
ues of M co id derived from CO observations for the galaxies 



in our samp le [NGC 3351: iDevereux et ail j 19921), NGC 3627 
Regan et al l (T2002J) . NGC 4501: IOnodera et alj<l2004h . NGC 4536 
Jogee et all j2005l), NGC 4569: iBoone et al ' 



20071) . 



NGC 4579: 

Garcfa-Burillo et al J J200gh l. Taking into account the different sizes 
of the regions probed by the H 2 2.12 /mi and CO observations, we 
can conclude that both measurements yield consistent results. 



5.3 Stellar and gas mass distributions and its implications for 
the M BH estimation 

Modelling the dynamics of stars in galactic nuclei is one of the 
preferred methods to probe the central potential of inactive galax- 
ies and directly determine the mass of the central BH. How- 
ever, in dynamical modelling, the total enclosed mass is con- 
strained, including not only stars and the putative central BH, 
but also other components that could be present in the nu- 
cleus, such as dark matter (DM) and gas. The estimation of 
the relative importance of these components is essential, since 
it will be reflected in the derived mass-to-light ratio (M/L ra- 
tio) and Mbh- The need for including a DM halo in the dynam- 
ical modellin g has been the focus of a number of recent publi 
cations (e.g., iGebhardt & Thomas! 120091: IShen & Gebhardtl 1201 ' 
McConnell et alj201 lbl:ISchulze & Gebhardtl201 lMGebhardteT 



201ll : |jardel et alj201 lilMcConnell et alj2012tlRusli et al.l2012al) 

It becomes an issue if the sphere of influence of the BH is not well 
resolved (Rusli et al. 2012a). The omission of a separate DM com- 
ponent forces the DM to be incorporated into the supposed 'stellar' 
mass. Hence, the stellar M/L ratio increases. Since M/L ratio is as- 
sumed to be radially constant this results in an overestimation of 
the central 'stellar' mass that, in turn, is compensated by a lower 
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Figure 11. De-projected enclosed molecular and stellar masses measured inside a given radius for the galaxies in the sample. Blue open circles and red 
filled circles correspond to the molecular mass measured from the low-resolution and high-resolution data cubes, respectively. Green squares indicate the 
enclosed stellar mass distribution derived for NGC 4501, NGC 4536 and NGC 4569. The error bars reflect the 3o" uncertainties in the H2 emission-line flux 
measurements. The dotted line in each panel indicates the Mbh of the galaxy derived from its sigma (Table[4) and the dashed line in the upper-right panel 
indicates the values of the Mbh of NGC 4501 derived from the dynamical modelling. 



Mbh- Significant amounts of gas could affect the measured Mbh 
similarly. If the gas adds mass outside the BH's sphere of influence 
its omission could, in analogy to DM, force again a higher M/L ra- 
tio and low Mbh- If instead, the gas is concentrated near the BH, 
its mass could go directly in the M B h, biasing it too high. Here, we 
will address the question of how much the presence of gas could 
affect the dynamical Mbh measurements. 

In the previous section we showed that it is possible to esti- 
mate the mass of the molecular gas from the H2 2.12 yum emis- 
sion line flux. We used this result to derive the de-projected gas 
masses enclosed within a given radius for the galaxies in the sam- 
ple. For this, we assumed that the gas is contained in a thin disc 
located in the plane of the galaxy and measured the H2 emission- 
line flux in spectra extracted from the continuum subtracted cubes 
using elliptical apertures. The size and shape of the apertures were 
determined by the inclination angle and semi-major axis of the cor- 
responding galaxy (Table O, in a way that correspond to circular 
apertures on the plane of the disc/galaxy. Fig. QT] shows the de- 
projected enclosed masses as a function of the radius measured for 
the galaxies in the sample. The molecular gas masses measured in- 
side the smallest aperture, r = 0.2" (or 10-16 pc), vary between 
4 x 10 5 - 8 x 10 7 M©, and reach ~ 10 s M© inside larger aper- 
tures. Particularly noticeable is the relatively high gas mass content 
of the galaxy NGC 4579, which shows the highest measured mass 
value of the whole sample inside a region of just a few arcsec (i.e., 
r = 2"). Although this value is slightly higher than the one derived 
for the 3" x 3" region of this galaxy (see Table [3}, these measure- 
ments are consistent within the errors. The difference is due to the 
different approaches followed in the determination of the contin- 
uum of the spectrum in each region. 

The ranges of gas masses found in the nucleus of the galaxies 



Table 4. Estimated black hole masses for the galaxies in the sample. 



Galaxy 


(o- 2 + V 2 ) 1 ' 2 


Mbh* 




[kms- 1 ] 


[10 6 Mq] 


NGC 3351 


100.9 ±9.1 


4.0 (2.6-6.0) 


NGC 3627 


122.4 ± 7.3 


9.8 (7.4-12.8) 


NGC 4501 


157.8 ±36.9 


31.5 (92.9-82.6) 


NGC 4536 


103.1 ± 14.5 


4.5 (2.2-8.2) 


NGC 4569 


114.2 ± 13.6 


7.2 (4.0-12.0) 


NGC 4579 


140.2 ± 16.0 


18.3 (10.5-30.1) 



* The numbers between parenthesis correspond to the lower and upper 
limits of the Mbh given by the errors in sigma. 

in our sample are comparable to the range of M B h usually found 
in the centre of galaxies. From the stellar velocity dispersion cr and 
the rotational velocity V profiles available in the literature, it is pos- 
sible to estimate the M B h of the galaxies using the M B h _ cr e (as in 
iGiiltekin et alll2009l) . Tableg]lists the values of tr e = (cr 2 + V 2 ) 1 ' 2 
derived for the galaxies in our sample a nd the corresponding Mbh 
assuming the M B h - cr relation given bv lMcConnell et al.l j20lT3> 
for spiral galaxies. The estimated Mbh are about an order of mag- 
nitude higher than the mass of the gas contained inside a region of 
r = 0.2" for all the galaxies except NGC 4579 (see Fig. [FT}. As 
we mentioned above, the latter shows a particularly high content of 
gas. The galaxy NGC 4501 harbours the largest BH mass from the 
all sample, which is consistent with the value M B h = 2.7 x 10 7 M© 
obtained from the dy namical mode l ling o f the stars following the 
method described by iNowak et alj 1 2007b, which is bas ed on the 
Schwarzschild superposition code of Thomas et alj|2004l) . The full 
analysis of the stellar dynamical modelling of NGC 4501 will be 
presented in Erwin et al. (in preparation). 
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But how large is the gas mass in comparison to the stellar 
mass, and how much can it affect the M B h dynamical estimates? To 
answer this question we have included in Fig. QT] the de-projected 
spherical stellar masses as a function of the radius for the galax- 
ies NGC 4501, NGC 4536, NGC 4569 and NGC 4579. These dis- 
tributions were obtained from the luminosity profile and the M/L 
ratio derived for each galaxy. The luminosity density profiles were 
obtained by deprojecting surfac e-brightness profil es using a modi- 
fied version of the algorithm of lMagorrianl dl999|). taking in to ac- 
count the effects of PSF convolution (see lRusli et alj|201ll) . The 
surface-brightness profiles (typically in near-IR bands such as H 
and K) were derived from a combination of ellipse fits to isophotes 
and 2D modeling of images, using high-resolution HST images for 
the smallest scales and SDSS or Spitzer IRAC images for larger 
scales; in cases where only optical HST images were available and 
strong near-nuclear dust lanes were present, we also used A*-band 
images derived from our SINFONI datacubes. Details for individ- 
ual galaxies, along with the dynamical modeling, will be discussed 
in Erwin et al. (in preparation). In the case of NGC 4501, the M/L 
ratio was obtained from stellar dynamical modelling. For the other 
galaxies (NGC 4536, NGC 4569 and NGC 4579, for which dy- 
namical modeling is not yet complete), we derived the M/L ra- 
tio frorriJhe_coloursof the galaxies and the trans formations given 
bv lZibettietai] ( 120091) and lBell & de Jond d200ll) . These M/L are 
based on scaled Salpeter and Chabrier IMFs for iBell & de Jong 
d200ll) and lZibetti et aljj2009T) . respectively. Recently, evidence for 



Thomas et al.l20 11 


; CaDDellari et alj2012l;IConrov & van Dokkuml 


20121; button et al 


|2012|). Hence, our M/L should be considered 



as lower limits for the actual stellar masses. From the analysis of 
the inner 5" of the galaxies in WFPC2 archive images, we derived 
average colours of B - V = 1.0 for NGC 4536, V - I = 1.2 for 
NGC 4569 and V - 1 = 1 .3 for NGC 4579. These values are in very 
good agreement with the ones report ed in the photometric compila- 
tion o flPrugniel & Heraudeat] i ll 9981) for the smallest (~ 6.6") aper- 
ture, and correspond to MfL K = 0.6, M/L H = 1.2 and M[L K = 1.0 
for NGC 4536, NGC 4569 and NGC 4579, respectively. 

It is clear from Fig. QT] that the gas content of these galaxies 
is negligible compared to the stellar mass, which exceeds the for- 
mer by more than an order of magnitude. In all the galaxies, the 
gas mass is no more than 0.1-3 per cent of the stellar mass, and 
even though the galaxies have considerable amounts of gas present 
in their nuclei, it is not enough to affect the dynamical Mbh mea- 
surements based only on the observed stellar features of the galaxy. 



6 SUMMARY AND CONCLUSIONS 

We have presented new AO-assisted AT-band SINFONI/VLT IFS 
observations of a sample of six nearby galaxies, NGC 3351, 
NGC 3627, NGC 4501, NGC 4536, NGC 4569 and NGC 4579. 
This is the first of two companion papers analysing the gas prop- 
erties of these galaxies, and is part of a series reporting the results 
of our SINFONI survey aimed at expanding the number of reliable 
BH masses in the nucleus of galaxies by means of dynamical mod- 
elling. 

The SINFONI spectra of the galaxies in the sample display 
several H2 emission lines and, in some cases, ionized hydrogen 
and helium emission-lines. This, together with the high-spatial res- 
olution provided by the observations, allowed us to study in de- 
tail the morphology and physical properties of the emission-line 
gas located in the innermost regions of the galaxies, and the star- 



forming regions associated with it. Additionally, we have derived 
the gas mass content of the galaxies and compared it with their stel- 
lar masses in order to check if it could affect the BH masses derived 
from stellar dynamical modelling. Our main results are summarized 
below. 

(i) A rich diversity of morphologies is seen in the H 2 emission- 
line gas distributions of galaxies, including bar- and ring-like struc- 
tures, centrally concentrated emission and off-centre irregular dis- 
tributions. When detected, the emission lines from ionized gas (i.e., 
Bry and He 1) showed also a variety of morphologies, not necessar- 
ily coincident with the H2 emission-line distribution. 

(ii) H 2 emission-line ratios point towards thermal processes 
as the main mechanism responsible for the excitation of the H 2 
emission-line gas, independent of the type of nuclear activity of the 
galaxies. However, non-thermal processes are not completely ruled 
out, and seem to have a relative stronger contribution in the regions 
closer to AGNs. This scenario is supported by the measured H2/Bry 
ratios, which are higher in the nuclear regions of the galaxies. 

(iii) The high-spatial resolution provided by the SINFONI data 
allowed us to resolve the circumnuclear rings traced by ionized gas 
in the central regions of NGC 335 1 and NGC 4536. The analysis 
of the Bry equivalent width indicates that the star formation in the 
ring of NGC 3351 and NGC 4536 started 6.9 and 6.5 Myr ago, re- 
spectively. Additionally, the He i/Bry emission-line ratio measured 
in NGC 4536 is approximately constant, suggesting that the star 
formation started at the same time along the entire ring. 

(iv) We found that a rescaling of the H 2 2.12 yum emission-line 
luminosity by a factor B 1200 gives a good estimate (within a 
factor of 2) of the total (cold) molecular gas mass. 

(v) The galaxies in the sample contain large quantities of molec- 
ular gas, with total masses in the range ~ 10 5 - 10 s Mq. Never the 
less, these masses are only a 3 per cent (at most) of the correspond- 
ing stellar masses. Hence, BH mass estimates based on the dynam- 
ical modelling of the stars should not be affected by the presence of 
gas. 
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APPENDIX A: 

In this section we specify the location of the spatial regions from 
where we extracted the ID spectra analysed in Section|4] The ID 
spectra were obtained by summing up the signal inside elliptical 
apertures or the entire FOV (e.g., ~ 3" x 3" for the HR data). Ta- 
ble lAl l lists the centre, sizes of the semi-major and -minor axis and 
PA. of the elliptical apertures. The PA. is defined as the angle be- 
tween the positive x-axis of the datacube and the semi-major axis 
of the region (measured anticlockwise). The spectrum of the re- 
gions called 'Ring' of NGC 4501 and NGC 4536 corresponds to 
the region between two ellipses. In these cases we report the mini- 
mum and maximum semi-major and -minor axis. For the region '5 
- Ring' of NGC 3351 we summed up the signal of the entire FOV 
outside the elliptical aperture specified in Table lAll 
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Table Al. Position and shape of the elliptical apertures used to extract the ID spectra analysed in Section|4] 



_ , _ semi-ma or semi-minor „ , 

Galaxy regions Centre . . P.A. 

axis axis 

NGC 3351 

1 - Nucleus (0",0") 0.75" 0.55" 90° 

2- Wcloud (1.05", -0.15") 0.3" 0.3" 0° 
5 -Ring (0",0") 1.15"- 1.15"- 0° 

NGC 3627 

1 -Nucleus (0",0") 1.0" 0.75" 90° 

3 -N cloud (-3", 5.6") 0.87" 0.87" 0° 

4-S cloud (2.1", -3.75") 0.87" 0.87" 0° 
NGC 4501 

1- Nucleus (0.2", 0") 0.85" 0.6" 45° 

2 - Ring (0.2", 0") 0.85" -2.25" 0.6" - 1.9" 45° 
NGC 4536 

1 - Nucleus (0",0") 1.05" 0.7" 135° 

2- Ring (0.25", -0.125") 1.875" - 5.375" 0.5" -1.25" 90° 
NGC 4579 

1 - Nucleus (-0.15", -0.1") 0.8" 0.7" 135° 

2 - NE cloud (0.5", 0.8") 0.75" 0.4" 150° 
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